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This report contains a description of a computer program 
for estimating the parameters of a mixt.ure of multivariate 
normal distributions with unknown frequencies, aeans, and 
covariances. The basic equations for the procedure are pre- 
sented for the first time here, with their derivation omitted. 
An example with the results of the computer printout is de- 
scribed for an artificially constructed mixture of three 
bivariate normal distributions. The method of using the 
program and the Fortran listing are detailed in this report. 


VII 


CONTENTS 


Identification... 2... 2. ee ee 
Purpose. . 2. 1. ee ee we ew we ww 
Restrictions . 2... + se ee ee oe 
Method . 2. 2 6 se ee ee eee ow 
Usage. . 2 wu ew we ee ee ew ew 
Example - Artificial Clusters. ... 

A, Listing of Input for Example. 

B, Computer Printout for Example 
Fortran Program Listing ...... 


PAGE 


ee ee) 


— So SE TT 


A COMPUTER PROGRAM FOR THE 
MAXIMUM LIKELIHOOD ANALYSIS OF TYPES 


I Identification 
A. E 
3. Written by John H. Wolfe, January 1963, Revised June, 1964, 
C. U.S, Naval Personnel Research Activity, San Diego, California. 
D. Coded entriely in FORTRAN II. 


II Purpose 

Given m scores on each of N individuals drawn from an unknown mix- 
ture of multivariate normal distributions, the program gives maximum- 
likelihood estimates of the means and covariances within each type, the 
relative frequency of each type, and the maximum likelihood of the sample 
for a given number of types. 


III Restrictions 
A, The program uses 14,554, words in the main program and 47,256 
words in COMMON. One input magnetic tape (Unit #3), one output magnetic 
tape (Unit #2) and one tape for temporary binary storage (Unit #4) are 
required, 
B. Restrictions oa Parameters 
Numbers of Variables <$ 
Number of ‘types < 6 


IV Method 
The program solves the maximum-likelihood equations by one of four 
alternative iteration schemes, depending on a control card option. 
Suppose that m measurements have been made on N individuals. Let 
Xj be the it Variable for individual k. Suppose the population from 
which the sample of N individuals is drawn is a mixture of r multivariate 
normal distributions. That is, the probability density, f(x), is given 


by 
r 
f(x) = f Ago, (x) , where 7 AZ = 1, A,r, and 


s=] s*} 


ij 1/2 -1/2 ij (x,-M3) oid (x,-M5) 


o 
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Here 5 is the relative proporticn of type s in population, |o,**| 
is the detorminant of the inverse of the covariance matrix for typ« 
3» and Mi is the mean of the ith variable for type s. 


Let us define 8, (x) = as (x) / f(x). SUBROUTINE DENSITY calculates 


a_, €, and 8, for each individual, Define the “generalized sample monents'' 


s? 


ps P 
{ui jap! as follows: 
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If a subscript is omitted or set to 0, the corresponding term on the 
right hand side of the equation is omitted, For example 


N 
po F ue 2 l E 
ijao ija oN ket Xi jak 8, (%) 


SUBROUTINE MOMENT computes a table of (yPs } with b<a<j<i and px<s, 
ijab --"- - 


FUNCTION U (II, IJ, IA, IB, IP, IS) looks up the correct jab from 


the two-dimensional table created by MOMENT, even when the inequalities 
on the indices are not satisfied, 


The maximum likelihood estimates of the parameters ure those values 
which maximize the function 


N 


r 
L=,f, log f(x) - or Ac DD. 


Setting the partial derivatives of the likelihood to zero results in 
the tollowing equations: 
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The moments (uy5) can also be differentiated as follows: 
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Where bus = 1 if pas 


0 otherwise 
The derivatives of the moments are computed by 


r] 
u 
FUNCTION DERV (II, IJ, IS, IA, IB, IP) # —t e 
a 6 
ab 


P « Po. Pp , ab 
Where Coa MP, 80 rp? and ®ab o, for afomd. 


DERV calls on three functions, ONE, EM, SIG defined as follows: 


ONE (IA, IB) ® Sob 


nm ab ps ps 
Em ( II, "J, fA, IS, IP) ® <= 
( ’ ’ D) , ) bei? ) (4 jp us ) 


SIG (II, IJ, IA, IB, IS, IP) = uy (oF Z Nene) + eee 
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The funczion EM used with s=0 gives the second term in ° “ij 
an P 
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: : ous, 
and SiG used with s*0 gives the second term in 4 . 
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The maximum likelihood equations (£5 = 0 can be solved in several 


ways. One iterative scheme is 
Newton-Raphson iteration, which solves the . inear equations 


i) s 
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for 46, ° 
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On the next iteration Ob s ab . sen, . 


SUBROUTINE NEWTON computes the vector 
s 


j 
2 8 P. 
A(IAT, JAT) © 2 £,,/? O° 


B(IAT) @ -f, and the matrix of coefficients 


SUBROUTINE MATINV, a standard SIIARE routine, solves the linear 
equations for age, and stores the result in the vector B (IAT). 


SUBROUTINE RAPIISON computes the values of the parameters for the 
next iteration First it determines if any of the increments are so large 
that rp + aa <0 or >1, if so, tha increment vector B is shortened until 


no Ady moves re mote than half the interval between do and a boundary point. 


d led 
Pp 
That is, if AA,<o, ddp>- zg, and if BA,?0, AA coy, Arter the aor 


are shortened, they are added to the oid estimates of the parameters 
to obtain new cstimates for the next iterations. 


The main rovtine also plays a role in shortening the increment 
vector. If the new likelihood is less than the previous one, or if 
the determinant of one of tae covariance matrices as determined by 
RAPHSON is negative, then the increment vector B is shortened to half 
its previous value and subrout:ns RAPHSON is entered again, 


Several altemative versions of (£55) can be written, 


At the maaimum-lixelihood points 


s s s 
= 1. 
f50 s Uo 1 = 0, hence Wy = 
s s $,,s s s s s 
5 . M. . = M. . = M. 
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Substituting for ue and us in Sh ve have 


Bay ae 


fay * Mis 


When METH * | on the contro) card, the subroutine NEWT.“ iter- 
atively solves the equations 
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When METH = 2, subroutine NEWTON solves the oriy'nal set cf equa- 
tions, horeafter referred to as (2845 = 0) 
When METI = 3, the equations used are 
$s 3 ’ $s 
f "l/u, (,f09) aj- a/uy 2 0, 
3 
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All three types of equations have the same solutions but their radii 


of convergence may differ. 


For METH = 2, 


the function { of, 
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The partial derivatives of fj 
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ji can be summarized in the following formula: 


rte (0,5 - 6; 5) 


are then easily written as 
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where o<j<itm and o<b<ac<m, 


The above formulas are used in METH © 2, 


er ee 
lvoo = “o 
s s 
fog 7 47 
s Ss _ s 5,8 
ayy," Mag > ay FN? 
are summarized by the formulas 
s s s 
ifay * Mig Sie 7 Om Sy gl eye's 


shy a eee Sig) Spsl Spo 85 08in 
a ap 


- (l= §4b/2) - Sh0 {1 - 6. Flo 


jo 


+ (te 8 1g)? 


Py 
ia jb 


If METH = 1, the functions 


y * ajey)}s and 


* (1 = 854) (6540) * 8 485)1 


bt 7b)? * 


S a2 sS Ss 
a3f, eae afi; 8 Mo 
Hence —— ee ae —, ———— . 
p Ny s Pp 
3 Fab ui, : *ab ( o) q Sab 


ab. 
Multiplying by ue and substituting, we have 
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Thus the equations for METH = 3 are readily calcvlated from the 
equations for METI] = 2, Only the matrix of coefficients has to bs 
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changed ,fsimply by subtracting fis P yg 
e uy 88a 


Instead of Newton-liaplison iteration, a methoc of successive sub- 
stitutions may be used for finding Me and O45: 


Since f£ * = 0 = ie - wens » the value of M® for the noxt iteration 
oi i oi i 
. s s,s 
is defined as Mi = uj /uy . 
: s s S$ ss s, $ $,,8 
a - = N - - . 
Since se 2 0 Vin Novy Tuy Hy (94; MMs) we can solve 


s : F : ss 8, 
for 955 after first substituting Miu for re: 


$s s,s 8,5 s s,s 
-, hy . ry , e 
o.. ® 44 j/Me MM. , where Mi are the new values = Wy/u, 


ij ij 
‘Ine new values of dy can be determined by Newton-Raphson iteration 
using only the equations foo 2 ke -l#Q0} 
Lbirferentiation of these cquations leads to the systen 
s 
r a Uo 


r ps - 

or £ u és 0 2 -1 
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When the control card ‘ILTH 2 0, subroutine NEWTON determines the increments 
associated with successive substitutions. 


Experience with the pro,ram scems to indicato that the NEWTON-Raphson iter- 
ation schemes have very small radii of convergence as compared with the successive 
substitution methods. ‘ueorctically, however, the Newton-Raphson iteration should 
converge itore rapidly once within its radius. ‘Therefore, provision has been made 
in the program tor running a fixed number of iterations with successive substi- 
tutions so as to get improved initial estimates and then switching to Newton-Raphson 
methods for tht remaining iterations, If the control sard option is -l, <2, or <3, 
then a certain number of iterations by successive substitution will be used before 
Newton-Raphson iteration by methods 1, 2, 3 respectively, The cantrol card number 
{DUMP specifies the number of preliminary iterations, 


The subroutine INITIAL determines the initial values of parameters preliminary 
to iteration. ‘The proper determination of initial values is crucial to the success- 
ful convergence of any itcration method, The initial valves determined by INITIAL 
are quite crude, and the rescarcher may wish to write his own version of INITIAL 
a ftor some expcrimcntution. The present version also allows the user to specify 
his own guesses of the initial values of the parameters cf vertuin types. 


The subroutine INITIAL takes a sample of 100 individuals and subjects them to 
a crude clustcring procedure. for cach individual, a count is made of the number 
of other individuals within a "box" two standard deviations or a side around it. 
That is, for individual k, the number of individuals j is counted such that 
| Xin? Xi S55 for all i = 1, 2,...0, m The individual with the highest count 


is made the centroid of the first cluster--that is, his scores rre the initial 
estimates of the moans for the first type. 


The individuals in the first cluster are crased from the sample and the pro- 
cedure is ropeated with the remaining individuals, 


The (nitial estimates of the A, are all equal to l/r, where r is the number 
of types. 


The initial estimatcs of the covariances are the same for cach typy and ere 
equal to the covariances computed from the samp)2 of 100 taken as a whole. 


After each iteration the subroutine RESULT prints the current estimates of 


the parameter for each type. At the ond of the last iteration, the program PLACE 
gives the probabilities of membership in cach type for each individual, 


J 


Oost Avaiic.oie Copy 


These probabilities are: 
P(individual k §& Type S) « de 8, (x,). 
A few words shouid be said about the indexing used within the program, 


First of all, the indices of the moments do not range from o to m and o co 
r as in our 2quations, but from 1 to me] and 1 to rel. Thus the value of 


U(i2z, 3, 1, 1, 1, 4) * te - The moments are conveniently calculated by 


i ; a : 
setting 2. =X: for each Xk 


1 


and z = 1.0. Similarly G(1) = 1.0 


and G(2) is the relative density for type 1, 8, (x,). 


The values for PERS (K) # Ae) 


cov(r, J, K) = “5 and COVIN (I, J, K) ® 1) , 


k-1 


also AV(I,K) = Mi ° 


The routine MOMENT collapses the (uP jab? into a two dimension:} 
array. 


The single index KL is uniquely related to p and s and the singic 
index IJ'&N is uniquely related to the indices i, j, a, and b. 


In general, suppose we have an array indexed as follows: 


There sare M indices. The first index varies from 1 to N. 
Each succeeding index varies from 1 to the preceding index. 


Let S(N, M) # nurber of elements in this srray, 
Then S(N,M) © ,£, S(I, M-2) 


and S(N,M) = WoT ® (Net) 


Let IX(1)2IX(2)a .-.z IX(M) be a sequence of indices for a particular 
element of the array. The one-dimensional index of the element is 


MoI 
Ke le Eh, (1X(M-161)-29d)/5 
M 
or Ke le 2) S(IX(MeIe1)-1,1) » 


These formulas are used by the Function yto look up values of vijab . 
V Usage 
A. Input (TAPE unit 3, BCD-card images) 


1. Title Card in columns 1-72, an) alphanumeric characters. 
Pae Control Card 


_COLS: NAME DEVINITION 
1-4 MX Number of variables 
5-12 NX Sample size 
13-16 IRM Number of Types Assumed. If IRM = 0, 6 analyses will 
be done assuming 1, 2, 3, 4, 5, 6 types. 
17-20 ITIERM Maximum Number of iterations, If blank, ITERM is set 
to SU. 
21-28 CONV Criterion of Convergence which all parameters must 


satisfy hetween successive iterations, If blank, 
CONV is sct to .0001. 


29-32 IRUN #1 if every iteration is printed, = 0 if ontv the last 
iteration 1s printed 

33-36 MFTH =-=0 1f successive substitutions 1s used +1, #2, 23, 
if various Newton-Raphson methods are used. 

37-40 IDUMP The numter of preliminary iterations by successive 


substitutions before Newton-Raphson iteration for 
M TH = -1, -2, or -3, 


3. Variable Format Card. This is an ordinary FORTRAN Variable 
Format Card according to which the Jata will be read, 
4. Data Deck 
N sets of one or more cards per individual, 
5 Initia) Estimates of Paramctors(optional) 
a stinmate control ca 
Cols 1-4 = K = TYPE # (1 through 6) 
Cols 5-8 = A Proportion of population of type K 
(all 4 digits assuned after decimal 
point) 
(b) Means for type K 


8 digits per mean, last 4 digits assumed after decimal point. 


(c) Standard deviations for type K (same format as (b)) 

(d) Correlation matrix for type K (1 row per card, same format 
as (b)) 

(e) Estimate control card for another type, etc. 


A blank card terminates the rcading of initial estimates. If 
no initial estimates are to be read, one blank card must be read, The 
sets of initial estimates do not have to be present for all types. For 
example, a set of initial estimates for type £ may be followed by a sect 
for type 3 followed by a blank card. If any initial parameters are read 


for a given type, all parameters must be read for that type. For 
example, initial estimates of means for type 3 must be followsd by 
initial estimates of standard deviations and correlations for vype 3. 


If no initial estimates of parameters are read in, the computer 
will generate its own by a clustering procedure. 


Muitipie runs may be made at one time by placing one batch of 
input “cards [A(1) to A(S)] followed by another batch. The last batch 
of data must be followed by two blank cards. 


B. Recommended typical usage. 
Under most conditio.s, the control card A-2 will be blank in all 


but the first 12 columns. The composition of the input will be: 
Title Card 
Centrol Card (Cols 1-12 only) 
Variable Format Card 
vata Cards 
3 blank cards. 


Another highly successful set of control parameters is METH = -1 or -2 
with IDUMP = 40 


C. Output (Tape No 2) . 
e natural logarithm of ae likelihood is printed for each 
N 
iteration: 421 108 £(x,) - N( t ag > 1). If the likelihood on an 


4 s=) 


iteration is not greater than that on all previous iterations, the 
computer prints out"Iteration--diverges". The parameters for each 
type are printed out on either the last iteration or on every iter- 
ation, depending on the value of IRUN, At the end of the last iter- 
ation, the probabilities of type membership for each individual, 

d s®s (x, ), are piesa, If sense switch 3 is down, the computer dumps 


ouc — monment 5 aut at? and the matrices (9,5) and (04?) following 

subroutine MOMENT, the matrix A of coefficients of che Newton-Raphson 

iteration and the vector 1-£)5) following subroutine NEWTON, and the 

matrices A! and the vector {4 a following the subroutine MATINV. 
D. Time estimates. 


Fo: METH = 0, time per iteration in seconds is 
Ts 1.44 N (mel) (m2) (rel) 

For METI. = 0, time per iteration in seconds is 
T = .U075 (mol) (me2) (me3) (me4) (rel) (re2)N 


+ .00002{ (mel) (me2)r) °N 


A typical run with \ = 225, r #1, 2,3. 4, 5, 6, m= 2 required 
50 minutes by METH = O on the CDC 1604. 


E. Suggestions for reducing runnin me. 


1, By intuition, cluster analysis of varisoles, factor analysis, 
or any other means, reduce the number of variables to those that are 
really important for discriminating types. 


2. Use any good prior information from theoretical hypotheses or 
any other classification-clustering procedures to provide initial esti- 
mates, 


3. If you have a large sample, say 10,000, don't run it all at 
once, Take a small sample, say 100 to 200 cases at random and run an 
analysis using a weak convergence criterion, sav .Ul, Get an idea of 
how many types there are and some initial estimates from the small 
Sémpie, Run the entire set of data using these initial estimates on 
specified numbers of types. The hypothesis Ne that there are r types 


can he tested against the alternative, Heed that there ere r-1 types by 
the x° test ‘ 
Sar -2 log va »® 2(L, - Leep 


with d.f. = (mel) (me2)/2 
where L. = likelihood printed for r types. Thus the 


x? test of the likelihood ratio on the small sample may give a good 
guess as to how many types there are, A final test on the entire sample 
could be performed specifying r, rel, and r-1 types on the runs, where 

t is the hypothesized number of types. 


vi Exampie e Artifical Clusters 


To test this program, an example was constructed consisting of three 
artifical clusters in two dimensions. The points in each cluster were 
generated by a pseudo-random normal deviate generator. The character- 
istics of each cluster are given in Table 1 below, The results of the 
computer printout are summarized in able 2. The points are plotted in 
Pigure 1. The 7S points in Cluster 1 are designated by triangles in 
Figure 1, the SO points in Cluster 2 are designated by squares, and the 
100 points in Cluster 3 are designated by circles, Drawn around some 
of the squares, circles, and triangles are larger squares, circles and 
triangles, The larger symbols give the classification assigned by the 
computer program with 3 specified types. If a point does not have two 
syabdols around it, it was correctly classified by the computer, It is 
evident from the figure that most points were correctly classified, and 
the computer's type< ‘ave clear cut boundaries whereas the actus) 
wausters overlap i'- somu dogree. 


TABLE 1 
Characteristics o% Artificial Clusters 


Cluster 1} | Cluster 2 Cluster 3 
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TABLE 2 
Charac:eristics of Types from Computer Program 
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TABLE 3 
Likelihoods and x? for Numbers of Types 


Number Natural Logaritha x? (with 6d. f.} 

of Typcs ot Likel}it..«! e2(L, - Vpo}) P 
1 - 380, 48930 
2 - 358. 96408 43.04924 tna x 1076 
3 -340, 36400 37,20136 .161 x 107° 
4 -334,83078 11,06644 .863 xX 107 
5 -325.63847 18, 38462 .$34 Xx 1072 
6 -318,02872 15.21950 .18s6 x 107! 


The data cards for input are listed in Section B and the output is given 
in Section C. 


A previous run with METH = 0, IRM = 0 (not given here) gave the likeli- 


hoods for 1 tu 6 types. These are presented in Table 3, along with associated 


x? values, The rosults indicate that the hypothesis that there are only two 


types can be rejected against the hypothesis that there are three types; but 
the hypothesis that there are three types cannot be rejected against the alter- 
Native that there are four types. 


In order to save space :n this report, only a run with IRM = 3 is given 
in Sections 38 and C, 


The method used was METIi © -1, which has 10 preliminary iteratious vy 
successive substitutions followed b, Newton-Raphson iteration. The previcus 
unpublished run with METH = 0 took 45 iterations to converge, while che one 
presented here took only :7. The difference can be attributed to the superior 
convergence rate of Newton-Raphson methods, iiowever, other computer runs with 
IRM © 0, METH = -1 and IDUMP = 10, 30, ot 10 sometimes failed to converge at 
all, once the Newton-Raphson procedure was started. if the likelihoods have 
not coverged to 0.1 by the successive substitutions, then Newton-Raphson 
iteration often fails. Thus the initial estimaces must be quite accurate if 
Newton-Kaphson iteration is to work. So far the various methods -1, <2, -3, 
appear to work equally weli. A'1 three converged in exactly 17 iterations in 
the present example. The run reported here took 13 ninutes on the COC 1604, 


SECTION VI 


A. isting of Input for Example 
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ie Pee 2 
ei} eth? 

19M I1L7 
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044-128 
20820129 
oe 7049 
1'3etiua 
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Nr ae1i4 
0*8-160 
~164~146 
~0hae148 
082-008 
=079-200 
10M—t 36 
21.46-085 
otter oiba 
141-194 
0482413 
1c2°129 
253-164 
20%6+055 
o0Anall® 
3040-162 
e ON Fo 06 
20700144 
2742077 
049-197 
O4n0158 
914-436 
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227143 
hho par 
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47? 929 
35° 978 
77° 18A 
720427 
<)97 245 
2352910 
-219 3245 
55° 157 

ew ARs 147 
-)79 %ai 
\eryi 194 
i: 
Tu 32° 
71 94845 
wh 24) 
HAWWG 
14 337 
#115 2392 
“114 215 
-179% 331 
25 312 
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41°98 
35 955 
=103 280 
7292-038 
9? 1é4 
7295 (53 
-186 108 
18 084 
"215 128 
45 400 
49-""" 
-161 288 
7 143 
“111 982 
116211 
90 066 
34 112 
“113 133 
13196 
ii 69) 
“186 130 
106 240 
“103 16° 
40 178 
-)38 157 
11 224 
32 391 
“101 185 
63-027 
132 180 
62 052 
275 092 
107 042 
208 306 
38 108 
218 337 
113<150 
139 147 
27 110 
14 199 
60-056 
240 179 
161 2% 
141 084 
86 61 

7 129 
23@ 351 
172 200 
230 153 
160 025 
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98 

99 
100 
10) 
102 
103 
104 
105 
106 
107 
108 
109 
110 
11 
112 
113 
116 
115 
116 
117 
116 
119 
120 
121 
122 
123 
126 
128 
126 
127 
13a 
129 
130 
131 
132 
133 
134 
139 
136 
137 
136 
139 
140 
14} 
142 
143 
144 
148 
146 
167 


130939 
17 974 
214 107 
135 378 
224-762 
206 224 
133-012 

9-049 
56 056 
18io 325 
163 76% 
1AD 106 
261 163 
205 139 
193 126 
359 255 
118 187 
109 N74 
2.1 038 
23 083 
322 093 
202199 
$4-147 
tlio yi 
10-928 
1lz 040 
106 029 
2lt 093 
148 118 
23 ies 
6 027 
136 0863 
79 059 

2-020 
31-089 
139 992 
287? 217 
44-002 
261 185 
179 160 
96 079 
$8 009 
22) 186 
232 180 
222 117 
84-027 
146A 076 
135 314 
3:97 270 
25-095 


148 
149 
150 
15) 
152 
153 
156 
155 
156 
15? 
158 
159 
160 
161 
162 
163 
164 
165 
léu 
16? 
163 
169 
170 
171 
172 
173 
174% 
175 
176 
177 
178 
179 
180 
181 
182 
183 
186 
165 
166 
187 
168 
189 
190 
192 
192 
193 
196 
195 
196 
i9? 


203 

80 
175 
233 
245 
261 


23- 


1869 
218 
226 

96 
259 

31 
112 
365 


200 
110 
122 
225 
147 
141 
018 
022 
169 
261 
018 
197 
032 
14° 
246 


306-035 
242-010 


142 


106 


198-102 


107 


076 


24-0546 


146 
324 


0646 
046 


62-022 


163 


037 


$1-050 


162 


071 


175-163 


198 
199 
200 
201 
202 
203 
206 
205 
206 
207 
208 
209 
220 
211 
212 
213 
216 
215 
216 
217 
218 
219 
220 
zal 
222 
223 
226 
225 
BLANKOS? 
BLANKOC 
BLANKOO. 


SECTION VI 


8. Computer Printout for Example 


© TTegalyya 
Li geRaty os 
1veR@ation 
ireRalion 
ITERATION 
LTeepation 
iteealion 
lTekaATION 
LTeealion 
ITERATION 
TteRaTION 
TteRaTIONn 
TT BRATION 
ITERATION 
(TeRaTion 
ifemalion 
iteration 
TteraTyon 


treeudeoor 
» LIweLtTwaAn 
Be, LPR PeOAS 
4, CimbtTaoe4 
Liwee Pwuon 
LTWeLTWOAN 
7, (iki LIA 
LEPL EROOR 
LIRELTWOAN 
LIMFL Twonn 
LIMELTROON 
19, LIMELT HOON 
11 PIVERGES 

13, (tebl Te00c 
14. t teeL wan) 
15. tIeeLTWOON 
1h. LIKE PeOAN 
17. LYMEL 1 WOOD 


Ah hnhahahaAhnhan 


Ahaha aswn 


Tytes 
TYPES 
TYPES 
TYFES 
TYPES 
TYPES 
TYFES 
TYPES 
TYPES 
TYPES 
TYPES 
TYPES 


TYPES 
TYPES 
TYPES 
TYPES 
TYPES 
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THIS 
Yuys 
Tels 
TmIs 
TwIs 
TwhS 
TTS 
Tas 
THIS 
Tus 
Twts 
Tats 


Tats 
Tels 
Tats 
Taps 
Tats 


Sa “PLé 
Sa-PLE 
SaMPLE 
SAMPLE 
SamMPLe 
SarPle 
SamPLe 
SaMPLe 
SamPLle 
Sa@PLE 
SamPLe 
SanPLEe 


SémPLe 
SamPce 
SamPce 
SanPLe 
SamPuiy 


0+ 495732926003 
©. $71059406063 
©. 36746558E003 
0. 304405638003 
©. S6ye4agaregs 
©, 3969033926003 
©. 4527414260983 
oe. 349 Q1014F O03 
@, $462825 76093 
©, 3447692SF¢03 
@. $463815686603 
oe. 3468752hE003 


@.94173960718003 
ww. 54 5A4823E 003 
©.340372 726003 
© .340N6407F OOS 
©. 3408640 QF ONS 


TYPES 


MAX TMi MeL ERELTROCD anaL YSIS Or 


ARTIFICIAL CLUSTERS weTRoney 2 yPReUtYivaRY TTretatluys 


SAMBELE SIZE 3 223 
wUMMER OF VARTABLES = @ 
NUMWER OF TYPES 2 x 


LPTERATION NUMPER 17 


LIMELTMUOR OF « TYPES IN TmtS Sar eLr 2 7.94 $O4n VEO 3 


rHsaracreseis7ic: Ch THe wHOLE SamPLe 


wBaN’, 
14h adn 
STANDARN CEvTATIONS 
1,46 1.47 
COPRELATIONS 
Leung +1849 
21849 1.9) 99 
CHAROCTERISTICS OF TyPE 1 7 


~ 


Tm PROPORTION Ob Tub PORPILATIUY FOP THIS TYPe 9 1484 


MEANS 
1.19 93 
STAAN Aw PeviaTIOns 
1.:4 29s 
24 


Wa Fe 10) 5 
Bost Avallans ver 


CORRELATIONS 
2.0nu0 «523% 
25231 1.4000 


CHARACTERISTICS OF TvPE 2 


THE PROPORTION OF THE POPILATION FROM THIS TYPE © 1346 


MEANS 
220 =1.92 
STANDARD DEVIATIONS 
1.28 530 
CORRELATIONS 
1.0000 +2462 
22462 1.9001 


CHARACTERISTICS OF TvPE § 


THE PROPORTION OF TE POPILATION FROM THIS TYPE © 1.79 


means 
1.41 1.79 
OTANDAKD DEvIaTiOns 
63 $.32 
CORRELAT ONS 
1.000 07498 
ty} 4.900y 


PROwS AT TTIES SF TYPO MEMAERSETP 


Be : AS 

\ . ot wurte erat 

¢ 0?) ge . 

S) 27490 wrt, ’ 

4 .v5? : 

5 , 29 vis : 

a « Ak .Y 

7 . 29 ee) ‘ 

“ 0 At ers . ‘ 

i) oe. mt ,wen i 
‘ pag Ree 
a 2A mes wade 
te 2 a8 a9 : 
+8 wore 275s . 
aq ohe A . 
"5 es eerie . 
ria) 29 2Ves ou 
v/ et 98S . 3? 
oa) 2 ro vcs ge ane 
1 7 2 a7? 
= 2 oh ais . ' 
2% a. Pt ,ved . 
ec . 34 9dr eu 
23 Femeee| 2e7t af Se 
aa go Ss 5966 F 
25 es ~ 9A? . 
ris . ts 92% . 
2) Sac aA ; 
an ot RY 4? eee 
24 P| af oa. 
3 oi $a Pela potas 
3) ony ovis os 
3é nF gat , 
as 0156 1ats as be) 
a4 1 7t 2979 ot! 2 
3? oA 0990 yo 
Jn fA? oa . ‘ 
v/ 179 OTs ‘ 
Ar) euyt 699% ace hat 
39 1496 2».4 o yf 
4, esd vk ell 

i eA .974 ede 
4? iit ~4Ae we RS: 
a3 ous 95% ett 
aa er da 1 146 as 
as of Ba ae ie eo! 
a4 0 h4 o9tn eo 
a) et 3? dA! . of 
aan 07 18 7A » i 
au «38 2946 . 14 
5, . #7 934 . 
S41 .19 9A ood 
§2 Pea | Pe el eo 
& 4 wtaa 94" ’ 
$4 «, 7n 2970 oa 
aA ot? VAL ed 
Ra 69 ot . onl 
a? ott? Woe eed 
54 015? Ab? . 
69 7799 Je . 1 
é e844 o ORS . 


26 


CJ 
© 
@) 
or 
Be 
< 
5. 
ay 
Ss 
CB 
) 
© 
Oo 


7 44 
»uné 
eU3A 
2349 
017 
ould 
213? 
~74h 
0 BY 
0217 
2nn? 
0 G14 
2un4 
0144 
e714 
018A 
ohn 
a QUA 
+999 
0949 
«Gul 
22% 
O17 
0 GAN 
eonn 
©0145 
2o8e 
927 
«009 
1094 
e004 
904 
2001 
oOfin 
008 
2184 
oN? 
768A 
2963 
200A 
008 
o CLA 
2014 
0004 
196n 
20n? 
o0n4 
2918 
«CON 
Bin 
2168 
0070 
2479 
9958 
urL) 
394 
9 904 
0004 
2984 
090 
216f 
007 
eiGA 
e000 
eOit 
28nn 


“958 write 


{99H 2300 
964 -onu 
626 ar Bl 
.9Aa3 2n0u 
2984 e.My 
463 ofall 
«254 ore 
949 Party 
783 Pee ( 
993 0 1G 
986 e>OU 
~ 396 eOU 
855 20g 
984 ote 
0099 AS 
2Ony 98H 
ed «994 
2383 2 08 
ent -958 
7us8 1961 


04) L.a00 
09093 1.400 


und 1968 
00935 e911 
6009 oars . 
2000 0991 


0090 6004 
0099 .094 
2066 028 
200 099 
2900 1.09600 
2000 097 
O04 0044 


2000 098 
6000 91? 
oONt »nde 
0000 4.900 
2031 0041 


6900 084 
«0C0 6966 
a) -094 
«000 0040 


2609 998 
2000 009 
6073 +908 
260 1.000 


«600 +190 
«000 04313 


e00n Pay 2 | 

2000 «$23 

2999 0905 ‘ 
2000 064 

2000 1406 

2000 2009 

2000 096 

2900 “4046 

»009 084 

2009 A40 

2009 003 

2009 oA94 

6000 ie0nd 

0909 .989 

e200 noo i 


SR MAE Re ae 


1a? 


Jenn 
«399 
1.598 
99? 
Leas 
99? 
Le. 7 
oot 
te tin 
7 99R 
a 
477 
Sean 
Lei en 
Leta 
2999 
RA 
1.°09 
Len 
» ra) 
~¥8? 
1994 
0997 
Leena 
Leu 
2148 
Letra 
BGT 
2275 
9990 
o¥Aa 
«999 
> re) 
Tein 
Deluna 
Lewin 
Letun 
» ) 
1.004 
0999 
997 
Let OA 
oGA4 
ob 3A 
le tn 
oH NS 
2997 
6 99R 
LeuNnna 
1.5.99 
e9An 
2994 
Dougan 
0999 
o“6R 
Aa 
Le. -A 
ie’ vy 
994 
Deon 
Lernin 
5 
«978 
Lerrn 
Leena 
dena 


re os 
weow: 
_o™ 
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Best Available Cop’ 


215 
216 


«784 
Zeu0n 
Zeunn 

2154. 
Lenn 
1-009 
1-¢9A" 
1-60N 
Leuan 
Lela 

AVA 

6979 
LTeifin 
1.9R6 

2988 
1.vUA 

998 
Tein 
1.0" 

135 

2709 
1eCONA 

0044 
1.0Nn94 

o 504 

990 

098% 

833 

906 

353A 


1.000 


on? 


214 
50) 
eu 
on) 
wed 
03S) 
2JC) 
»o4) 
Pree 
«005 
0123 
£021 
0090 
0094 
245 
ouGs 
005 
0000 
2009 
865 
2291 
094 
959 
900 
494 
Ont 
948 
2167 
005 
462 
#00) 
903 


oil 
en 
«thd 
0200 
or An 
2700 
eh 
oe ANG 
27ng 
oe Ihe 
-n00 
“Aut 
900 
eng 
era 
e806 
600 
rod 
enn 
e100 
ead 
160 
2206 
enn 
er nts 
o)H0 
eff 


o00. 


100 
«900 
noo 
eng 
700 


29 


a ad 


yonteple CODY 


pra FS 


*% 


Yy CAls 


VII Fortran Progras Listing 


aan 


Re Se 


PROGRAM Tyee : TYPR0000 
DIMENSICR v¢S, 8 20766), PEREZ PERS(7) AVIS 07) .COV65,5575 oCOVINGSTYPEDN0R 
105 o 7) e MLPHAC7)  DETERMET), GO7) THES) o 1D64)0V61265 28) 0 AC126,120) BCI TYPEO0O? 
226) HVIT2H9 MECORD C22) 0SN05)  KGAM C5 5101) MO CLune dd TYPEQUOS 
DIMENSION emtei2) TYPE0004 
COMMON X67, PER, PERS AV sCOVeCOVIN, ALPHA, VETENM, Go Vo dot, AV o RECORD, SOTVPRO0US 
LoMY, MNT, TA, WIR, LRT LUMNE LAL NY DUNK» JATSM OC CONV, TTERM OL Up KEAM, MATL ETYPE0006 


CPue IMS TL eo TITER, PROB,METH, TOUMP TYP&0N07 
EQUIVALENCE (XSAN, VD, CAS MAT OND V¥PE9008 

1 READ INPUT TAPE 3.1%, CMECORNE Ts Leind 2) MUON, TRAM, | VERMSCONV, TRUNTYPEO DUD 
1.MOTM, IDUMS TYPB1020 
101 FORMAT C12AAs14, 18,214. FO Ae S14) TvPG0021 
1Fcwn) 23,9302 TY®E1012 
MYSNUMBER AF VARIAGLES. TyPko0lS 
TWE PROGRAM READS BATCHES OF DaTA UNTIL 2 BLANK CARLS ARE ENCOUNTETYPED014 
1860, TYPG0025 

2 SFC TERMS, 3,4 TvPEv026 
ITE@MEMAXTMUM NUMMER OF PERMISGABGLE ITERATIONS, TYPE0017 

3 TERMS: TyP&0018 
4 1F¢CUNVIS. 8.6 TvPR0029 
CONVOCRITERTON OF CONVERGENCE uMICH ALL PARAMETERS MUST SaTtSty eeTYPROORO 
1TWEEN SUCCRSSIVE ITERATIONS. TvPG00as 

3 CONV® 0.001 TvPkO0Rs 
6 MNL)8MKel TvPEocENs 
REWIND 4 TvPe0eR4 
XNKANK . Tvegooas 
NEVORINTECWNE/10U. 01 TyvPéo0z6 
LYMNXOCMNLAEMNLOLINCMN LEZ ACMI OS) II2G YvPRoce? 
Kannec svegoesé 
Lyen TyvPeo0s® 
READ INPUT TAPE 3,152.06 uTC1), 189512) TVPE0030 
0? Fommatci2aas TvP960033 
7 CALL OLOCH ENE, 3009, KAMN, KET, REND ,LONG) TvPseess 
00 @ Jel,LAne TvPeeess 

6 READ INPUT TAPE SFMT, CHET odd TOL MHD TVPB0O34 
DO O Jel Lara. nev TvPRCOSS 
Lee hes : TvPeeese 

DO @ leisMy TvPBe037 

@ RGAM CLoLuvemiteu? TVPROC36 
RRAmeSanei@ CF UP TO 1uC CASEN USED FOR INITIAL ESTIMATES OF PaRanTvPEENS® 
1ETERS. bala (LJ 
WRITE TAPE Gol ORCL dds lat eMEd.. 182 +L ONE) Tveenes 

16 cme, i367 tvPeeeds 

du Cary Enittar TvPGeees 
Tastan thle Li bas 
TRMONUMBER OF TYPES ASSUMED. 17 Immed, © DIFFERENT ANALYSES ane TvPEeees 
DONE ASSUMING 1 TO 6 TYPES, TYvPRES 46 

ie ¢hamyi1,13,592 TvPH0047 

41 aaphes TvP60048 
a2 [Artes TVOECECS 
meTwenatn rvegeese 


XTReIR PYPE 0053 


RL KeCIM1ec1R109e1)I972 TYPE 0052 
VAT © (Mer eCexrezr #1@I79 TYPE OD ySS 
PERS(K)SPRAPORTION OF POPULATION OF TYPE &. TVPG cCS4 

DO 933 wer,181 TVPEC OSS 

16 €12m993.7.52922 TVPE 9056 

O42 16 (PERG CKI113,13.023 TYPE 0057 
q3 PERSCK) EL. /XKIR T¥PE HOSE 
073 CONTINUE T¥OE 0059 
ITER|N TYPE 9060 
Ipives TVPG 9061 
PROBAR01N4A5 76, 008 TVPE 0068 

24 1TEROT TERS, TYPE 0063 
IE CDETERMI12d SSyro3Sn163S0L TYPE 9064 

3$00 DETERMI 1)Ien.0 TVP6 9963 
GO YO 346 TYPE nee 

S901 CALL MOMENY VVO6 0067 
WRITE CUTRHT TAPE 2.9916, ITER, 18, PROD TYPE 9 00e 

0016 FORMATE 10m ITERATION 13,15. LINELIMAOD OF 13.23 TYPES Sy THis SATVPE 006 
1MOLE S61A,a) TVPG 007¢ 

16 (PROBPRABA) 3406, 547,347 VYPE 9073 

346 IDI velDive, TVPB 0078 
LTERASTTER.: TYPE O07. 
WRITE OCUTPLIT TAPE 2.9915, 1 7ERA TVPG 0074 

9025 FoRmavcifw TTEMATION 13,0m DIVERGES) TvOE C079 
DO 2445 war, yatmn TYPE 0076 

3445 BCKde =|, SeRtK? VV¥0G0077 
16.101 v—099941523646,315 TVPG0G 78 

3446 DO 3447 Rae, Jaron TV¥HE 0078 
35447 BEKre eBtK) TVPG 0080 
Go TO 218 7V¥9E 0001 

347 IDI veo T¥PGe082 
16 CmETMD 304, 340,308 7¥960 083 

3446 TECTTERST DIMM) 340,348, 9406 7¥06 0004 
3444 METUeemETH TvPE 9009 
TpuMPs THE NUMMER OF PRELIMINARY LTGMATIONS By SUCCRESI VE Tv900086 
SUGRTITUTIANS REFARE NEWTONCRAPHSON ITERATION FOR NEPNO>s ,-8,0R oe ZTIVPR NOS? 

F440 PROQASPROA Vv0Gg00e 
IFC SENSE SuITCe 39540, 2469 TYPE 0088 

240 URTTE OUTOT TAPE 2.9910 TYPE 0080 
90g" FORMATE Ze BQNENT) TVPG 0091 
WRITE OUTHHT TAPE 2.0014, CEVEL Sd o ISL ONaMds 104, 5 MUN D 71069008 


WRITE AUTHNT TARE 2,090146, COCCOVE Le rod, 00H) op SOlo ME) oOo TRG) TYPE OSS 
PRIVEE CUTHHT TAPE 2:0OL SG, CECCOVING Eo toM do LOL oN de dOl MND oROge IA) TYRE og04 


240 CONTINUE TVPG 0009 
194 Care NSafOn TVPE OCDE 
IP CQENSE SulTCu3) 399,49. TVPG 9087 

Wn URITE HUTA? TAPE 2.0012 TVPG 0000 
90q1 FORMATI7™ SBUTOND TYPE 0009 


WRITE NUTHLT TAPE 2,99016, CAAT Jd of Sis JATOND VEL, UATAN) TYPBQL0C 


32 


WRITE OUTONT TAPE 2,9914, 0801). 102, caTMn? 
FORMAT *6F 49.6) 

CONTINUE 

TFCMETH) 47155115,948 

IF mETH @ pn, SUCCESSIVE SHASTITINTION as USEN. OTHERWISE 
NEWTON*RAPUSON ITERATION, 

Catt MATINGCA,JATMN,O,20NA) 

TF CSENSE SwiTCw3) 551,284 

WRITE OUTPIT TaPE 2,9912 

FORMATOC7M wmaTINV?D 

WRITE OUTPUT TAPE 2.9914, COATT Jo LO po JATMND JO, UATME) 
WRITE OCUTPIT TAPE 2.9924, (BCID, Ter. Jaen? 
CONTINUE 

Cate RAPHSAN 

TF CTRUND 19,16,18 

TF TRUNGO.ANLY THE LAST TTERATI‘N 18 PRINTED. 
Cart RESUL? 

TE CTPERTTBRMI17. 19619 

DO 26 med, jaTey 

TA-O(K) 

TASAGSF CTA, 

IF CTASCONV) 10,16044 

CONTINUE 

SOC TRUND 29,2022 

Cate RESULY 

Cale PLace 

16¢38922.9901 

1F(1RO919,408 

Gnd FILE 2 

Catt Exit 

OND 6001000000) 


ie Lh 
a 

rvegeies 
TVPG oi10 

TYPE 9408 
TYPE 9406 
TVPB 0207 
TV¥PB C108 
TYPE 9209 
TvPB 0110 
VVPB 9213 
TVOR 0212 
AL. T3et) 
rwete128 
TVPG 0933 


cise 
198s 30 
rwestige 


aqgaaa 


aaa 


aaaa 


SuBAGITINE INITIAL VVOE°LS2 
DIMENSTON 765, 39°3795766), DERI7)  PEREE7) cAVESs79.00V65,5.7) -COVINGSTYPE NSS 
1065e7),ALPMAC7TI,DEVTERMI 7D, GO? IKE IDC 4b oVE1265 28) 0 AC126, 120) BCL VPEGLS4 
226), BV61 2H) MEFNRNC IZ) SACSDKCAMESe On mMaTeL nelnid TVPE 0335 
COMMON 017, PER. PERS AVeCAVACOVIN, ALPMA, DETERMsG.V2 ho Be RVs RECORD» SDI VPE 436 
LOMM MRL IR MIRe T Roo] UMNK KLM ONT e MNKs JATMK CONV, | PERM LK RQAMe MAT ETYPR NG 37 
2Pm,IXSID,TTER.PROR. METH, 1 DUMP TVPE 0338 
EQUIVALENCE tXSaM, Vv), CA MATSN) TYPE 9139 
TIS SUBROUTINE DETERMINES INITIAL ESTIMATES OF THE MEANS AND COVATYPE D9 46 
TRIANCES OF TWE TYDES AY APPLYING A CLUSTEMING PROCEDURE Th A SAMPLTVPE 99 44 


2€ OF UP TO +n: CASES. TYPE9242 
FIRST THE SAMPLE MEANS AND COVARIANCES ARe DETERMINED, TYPEQL43 
MUxekx TVPE 0344 
DO 22 1e1,4x V¥PE43 45 
Av(teties. TYPE C946 
BO 1 Kei.iv TVPB 0947 

1 AVOTOLIGAVE Te TP ORSAMET OK) TVPG 0148 
22 AVE Ted ebvelss ds X TVPE0149 
DO 3 lei.My TVPE0LS0 
0O 3 JelomMy TYPE 04S3 
COVE odetda’. Add TERT} 
DO 2 Ketoly VvPE0LS3 

2 COVED ode LMCOVEL Vet doKRSamel KrerSamcy,K) TVPB0254 
COVE Lo del PMCOVE Le uet PML we AVC LT, LICAVE NSS? TYPE0SSS 

3 COvedel,2eCOVEt.J.1? TVPG 0156 
DO ¢ le1.my TYPE 0387 


AN LK BY Ly MATRIx 1S COMPLTEN, AN ELEMENT CORRESPONDING 7O A PAIRTYPE0LSO 
+ OF POINTS 18 + IF AND ONLY IF BOTH POINTS LIG WITHIN & Bax wNQSE TYPE 9989 


2 SIDES AME GONE STAs Marr neo 'arres LONG, TVPB 0860 
SOveCOvel.1.1) TVOG0163 

@ SDCr»sgoRTECSDV: TvPB 0162 
DO 7? Keqety TVPG 0363 
00 ? Lemoty TVPB 0904 
DO S Lei .my vvPEG 0965 
RaNSAMEL My oMSAMCT OL? TVR 0166 
ReARSF CAR) TVPG 0907 
16 (MeSN11)9)52656 TvPR0100 

6 CONTINUE TVPE0269 
MaTCL on) 84 TVPE 0170 
MATEReLIOY TvPG0871 
GO TU 7? TvPE 0872 

6 MATL oRDON Tv0G0979 
WaTemeLder TvPG 0174 

7? CONTINUE TVPG 0878 
DO 34 ™s?,7 TVPG 0176 
THE CENTROTO OF A CLUSTER IS THE POINT wife Tee QREAVEST yUNBER OF T¥PH 0877 
TOMER POINTS 18 & AON ARAUND IT. TvPG 0170 
FIRGT THE INDEW OF TwE PRINT OF GREATEST VENSITY IS FOUND, THE TvP TYPE 0179 
16 MEANS ARE TRE CH-ORDIAATES OF 14E DENSEST POINT, TvPG 0186 
Mana TVPG 0183 
VO 1! elo TvPGosea 
L@en TvPG0863 
DO A Letely TvPE 0104 


aan 


aac 


= 


= 
= 


1? 
i) 


14 


a8 


4? 
122 


LL) 


“93 


48 


a0 
te 


iRT CLUSTER TmUs ILI 


LBat Romatcy ,«) 


TYPE 9202 


LE Cl MeManrys eo’ 69 TYPE 0186 
MaxelA TYPE 090? 
LAmen TvPE ogee 
CONTINUE Tv¥PENL8S 
DO 11 Iei.uy TYPE 019C 
AVC To MPCIS Amel LAK) TV¥PE 9993 
DN 31 Jelouy TVPE 00S 
COVE toto MPaCOVE Soto 9? TY¥PE 193 
COVE Lo deMI COVED TM) TVPE 0394 
DO 34 was, i TYPE 9908 
TE COAT Ct 09914, 40 92 TYPE 9996 
DO 13 Lai, TYPE 997 
MATERoL De, TYPE E300 


ALL PREVIQUSLY CLUSTERED POINTE ARE ERASED. THE CENTROTD OF THe NETYPEC O00 


CONT t hue 

0O 1 fet,ux 

DO 13 Jei.unx 

Actevie Chvtletoi) 

CALL MATinu( A, Mn .R.0.DA) 
ADsaRSFcHarsna 

DaecSOatrh¢aps0Ay 0A 

Do 16 «(Ke4,7 

PERSUK D004, 

DETERMCKIS DA 

DO 1@ Jag ,mn 

06 1% yer.mm 

COVINGE dom deacteud 

THE ROUTINE ALSO READS IN INITIAL ESTIMATES UNTIL A BLANK CARD 
1¢ BNCOUNTRRED. 

RBaW INPUT PAPE 3,399,K.H8 

FommATCL 4,646.49 

16a) 10529.18 

Rees 

PERQIN EDA 

READ INPUT TAPE 3,285, CAEL ORD, Jer, mED 
READ INPUT TAPE 3,123, CONCT I fadremnt 
Fonwmar (SFA, 4) 

BO 1% yet, mn 

READ INPUT TAPE FZ.1NT.CCAVEL I KPe LOL. mNd 
00 *@ lei,ux 

COVES DoMeCOVEs  yeSdeSDelsesDedd 
ACL UdOCOVET ed K) 

Care Malinved. Oy. Bee DAD 
ADeahSFcDarsha 

Baw SOntl ¢cans0Ar0ad 

OETEPMEKIeNA 

BO #0 tei.en 

BO 20 Jei.me 

COVINGL ed erehetoug? 

@0 70 37 

RETURN 

ENDG ode Dene? 


AE THe POINT wht fog GREATEST NuMOgR OF PREVIVPE Rca 
ZIOUSLY UNC: USTERED POINTS wITHIN A BOX ARQUND IT. 


TVPG 0203 
TVPEeR08 
TVPG e803 
TvPEAgC4 
TVPE CROs 
Ute Lt 4 
TvPG ese 

TVPR R98 
TvPR eg? 
eit 

e 
TvPsedae 


SUBPOUTINE DENSITY TYPE 9240 
DIMENSION WS, 35077, 7209), CF RO7) e PERSI 7) cAVE 507), COVIS, 507) -COVINCSTYPE D244 
DeSoPVs MSc Pua 7) NETEAME FI, GE, XCD e ID CHD oVEL2N, 28) 0A11 26.126), BIL TYPEQ 242 
226) BVE1LZ6Y MEGNAN Dre SNCSIAXSAMESr gods MalCingerHed TYPE 9243 
COMMUN X27, PER. FERS, AV» COVeCOVIN, ALPHA, DETERMLG, Ve Ae BsRVe RECORD: SDIYPE D244 
LeMXMXZ TR MIRO TR ye LUMNE KL X eNO XNX Ss JATMX CONV, LPERMAL NO XSAMSMAT LET YPEQ 245 


2PHL IX, ID? [TEA OROG, METH, IHUMO TYPE 9246 
EQUIVALENCE (XSAM, VD, ¢a, MATIX) TYPE D247 

C FOR CURRENT ESTIMATES OF TWE PaRAMETERSs THIS ROUTINE COMPUTES. TYPE 0248 
c ALPHACHM MMOL TIVARIATE NORMAL DENSITY FOR TYPE K AT POINT 7, TYPE 0249 
Cc EPHaCENSI Tv AT 72 FOR THE MIXTURE OF DISTRIBUTIONS TYPE 230 
c GOK) SALE MAsEPH AT POIAT 7. TYPE 0251 
Do 3 Ke2, tai TYPED 252 

ALan, TYPE0253 

DO 2 Lei.My TYPE 9254 
Lasley TYPE 9255 

DO 1 JajoMy TYPE Q256 

L ALBAL@COVIN CUT Ma eC 2 OT are aVET KPI REP COL wmAV EEK)? TYPE 9257 

DP ALBALON SHOVING Le Lem re reT aro ave lM) ) med TYPEQ258 

3 ALPHA (K)2 DETERM (KIKEXPF CAL) TYPEQ 259 
EPnWed. TYPEQ26U 

DO 4 Ka?, ity TYPE 926% 

@ EPHREPNOPERSIK) wal UMA(K) TYPEQ262 

DC 8 Ke2.Tay TYPE QZOS 

S GUK@ALPHAIKI/EDH TYPE 0264 
Gei)e1. TYPE 9265 

Cc THE ADDITIAN OF A DUMMY TYPE OF CENSITY Gay SIMPLIFIRS THE COMPUTATYPEQ 266 
Cc 1TION OF LOWER ARCER MOMENTS WITH ONE OR MORE TYPES OMITTED, TYPEQ0267 
RETURN TYPED 268 

END Lede Cer. 1? TYPE 6269 
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Best Availasle Copy 


rad 


aaaa 
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SUBROUTINE MOMENT TYPE0270 
DIMENSION wt5.3°5720706),PERC7) PERSE) oAV(5s7) oCOV65.507) oCOVINCSTYPE 0272 
1055 7s AL Punt 7d DETERMC7),GO79,1K04)01D6 4) 010126. 283 506126,226) BCL YPE C272 
226) HV 01269, RECORD(12) 0 SACS) NCAMESA1000 MATER, LUD TYPE0273 
COMMON Ko7, PER.PERS. AVoCAVeCOVIN, ALi*td, DETERM. GeV 2 hr 8, aVoRECORD. SDI VPEOR74 
LeMM ML TA KIRD TRL o LMNE KEL MO NM OMNK, UAT Rs SONV, FTERM LL RGAMS MAT ETYFEC RIS 
2Pm, Ik IDs TTER. PROB. MEH, [DUMP TYPE0R76 
GQUIVALENCE (2SAM, Vv), CA. 4ATOY) YVP&9277 
THIS ROUTINE COMPUTES THE GENERALIZED MOMENTS VETUMND KL) FOR VaR ATYPE0R76 
VOLES IToJom, AND N. AND TVPES K AND L. THE MOMENTS REALLY WAVE Six TYPE0R79 
2ENDICES BUT ARE STORED AS A PeHIMENSIONAL MATRIK, TREANGULAR INDEXTYPED200 


SING [$8 USER TO ELIMINATE DUPLICATION, : TvPEORO1 
DO 1 TumNs7, LUMNK TvPS0LO2 
DO 1 Kiel, wiX TyPECgES 
VET IMM OKL Oa. TYPE08O4 
Kanned TvPgoses 
Int*taclze TvPG0266 
REWtaD 4 TvPG0897 
PROBOXNE TvPR02E0 
00 10 Ke2,1@1 rvPgosed 
PRODAPROB-yNNaPERS CK) TYPEUROO 
CALL BLOCK NM, $000, RAMN, RST. REND, LONA) TvPReROs 
READ TAPE 4, CCMCT Ido TOL MNP oe e LONG) Ala ttyl} 
TAPE 4 18 MEAD IN BLOCKS OF 3070 Cases, TvPR 0803 
00 17 KAe1,L ONG TVPG E804 
2¢1904.0 TvPGO8OS 
THE ADDITIAN OF A DUMMY VARIAO Fel SIMPLIFIES THE COMPUTATION OF hit tb 
1OWER ORDER MOMENTS OMITTING OnE OR MORE VARIAOLES. tv9G0807 
00 3 leisy TvPEaROS 
ZeLotdenel Ad TYPER8O9 
CaLL DENSITY TyvPG0Se0 
PROSOLIMELIMOOD OF SaNPLE YVPEO301 
16(GPH)7,7,8 TvPEO308 
INDE MOKGTonAes TvPEOS0S 
WRITE QUTMIT TAPE 254200 INDEX, EPH TvPG03e4 
FoMMATIASH NEGATIVE PROBAGILI Ty DENSITY FOR COSERVATION 39,3 © TYPEEZES 
ad ae TvPGEse6 
eo? Ty¥0G0307 
PrupePRODoLOer tern) TyvPGese8 
eLen TvPaesC® 
IF CMBTHD 19,92081 TvPGogl10 
00 4 Ke2,1@i TvPhOS13 
bO 4 LeioK YvPG0318 
MLenLes TvPGesa3 
OT1eacmpearr) tve60314 
Tumned TvPR0919 
06 4 leiemyi Libe- be 
OTzeGTise2¢ 1) “ 4 He 
00 4 Jedel ; 

OT3eGtaezeu : ie 
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Un 4 Melos TrPe used 


Te 2GTsu7¢™) av EE See 
UN ¢ Nz’a™ Tyrtn 322 
LUMNal uM set TyPEUS23 
VOTINN ML I SVCT UM AL dearest} TyPenges 
CONTINUE TYPE S2> 
TB OKAMNIS SLC TYPEDSZ6 
DO ® Jumnsr.luvn< TYPEDS27 
DN s KLELetX% TyPeuds25 
VOL SON ML aves eet Z ONY TYPED 329 
RE TIAN TYPENSSO 
AL 2° TyPe 331 
DO 14 Kel.tnl TYPE: S32 
FRR Mb Yes: swe QE are fh aL Ter “MENTS TYPEDSSS 
LUMA 82 TYPEDSS4 
DO 14 Je2.¥! TYPENS3S 
GTie Ze}re wend TYPEDSS6 
JSes TYPENSS? 
09 14 Jal. TYPEUSSS 
VCLIMN KALI VETotNe Ye a Bara Wa eb) TYPESSSI 
JSzused TYPEY S40 
LUMNS TUM OS TYPEQ 341 
No 14 Laie TYFET S42 
VEL ALIBVET AL Vesa’ TYPEDS4S 
AL SKL ed TYPED S44 
GO tu 17 TYPEDS4S 
END ole ters * TYPENS46 
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FUNCTION Ue tds JAoMAONAS MAL A D TYPES 347 
DIMENSION wes, Piel do 20H), PEROT) PERSE7) oAV65 07) COVES, 507) COVINISTYPE D348 
LoS ede ALPHA sd, DETERMI70, GE79,1N64) o 1D64) eVE126578) 04(126,1260) BIL TVPEG S49 
C20), HVC126) RELIADECL 20. SNCS} MSAMESS ICCD MATTED. 419.9 TYPEC ISO 
COMMON MO7, PER. PEAS. AVS SAV eCOVIN, AL PHA, DETERM, GV eho Ve RECORH SDT YPEN3$S1 
LoMM MLTR MIR LAS LIMNE MLO NWO NK »pJATMM CONV, ILTERM LUO MCAMs MAT ETVPED3S2 


2PM. IDs TITER, PROB. METH, I DUME TvPE0SSS 
EQUIVALENCE (NSRM, V0, Cae MAT OND TYPEL 354 
THIS ROUTINE FINDS THE MAMENT FORRESPQNUING TO 6 INLICES aY COMPUTTYPEDYSS 
TING THE 2 TADICES DF THE MOMENT IN TWE MATRIN V, TYPE S56 
Iucqoata : TYPE03S7 
Incodega TYPED3S8@ 
Inc3o sma TYPENYS® 
Iecapena TYPE 360 
TwE FI@ST @ INDICES ARE PUT TN ORDER FROM SMALLEST TO Largest TYPED 3O1 
DO 3 Jele4 TYPED 302 
Mate10o TYPEDS6S 
OO 2 Ie1,4 TYPED 304 
TP CwATeITKCTDD Aedel TYPEN3OS 
Watelxci) TYPG0 366 
Kel TYPEC367 
GO tu 2 TYPEU3EO 
Cont INuE TYPE0 S69 
Ixcadege: TYPE1970 
Iptyrenat TYPE0S72 
TME LAST 2 INDICES ARE PUT Ih ARDER FROM LARGER TO SMALLER TvPkO372 
IF (Whol 474,505 T¥PG0373 
Mavala TYPE0S74 
LATena TvP?G0375 
60 TU 6 TYPE0376 
Katana TVPE0977 
Lavera TYPE0R76 
[umnet TYPE0S79 
THE APPROPATATE INDICES HF THE MOMENT MATRIX V ARG COMPUTED, T¥PE03AC 
DO 6 lei,4 TvPG02 32 
LORei TvPG03e2 
DO 7 Jeiel TvPRO SES 
LORACLPMACIDELPeQeudISU TvPE03O4 
LuMNe] JMNe) PR TV9E0SES 
RLecRaTecmatot sp s2eL ar TvPb03bO 
UeVELJMNe Ri) TVPE0387 
RETURN TvPR0S06 
En0(053.0,9,09 Ty¥PE0gS9 
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TYPE N3Oy 
TYPE 301 
TYPE 1392 
TYPE 393 
TYPE, 394 
TYPE 9395 
TYPE n396 
TYPE 397 


aa 


br NT Ine becitetu, tae bon 88) TYPE n398 
DIMESSIOS wes, t Ve7CAD, GEREMSCEQSI 7) DAVE Se 7) sCOVIS, Se 7d COVINESTYOE ~ 399 
Te Ge The ML Pal PME TEE TPL GEV, IKE MD SINGE So UC 1 Sh, ZBI AL 826,326), 803 TYPEN Gar 
PEA, HVE PHY REPRNANC FEL SMOCSD. GAMO SS 1 > 0,MOTEY 41? TYPE N49 
COMMUN KO 7 PER PERS AVEC IV ECIVING AL POA, DE TERM, Gr Veoh HV oRECORN» SNTYPE NGC? 
POM MKT LOUTH Lees Le Re EAM ATOM CONV, TTERM OLX MCame mal ETYPE 4,3 


PPuwr,T he I De TER, ORNR, METH, TA Me TYPE nana 
EQuUIVALES TE CKSAM Vd, Cho wato ed TYPE 4455 
EM TS a Team wel tw APPEARS IN Th Pel aL NEREVATI VE AE Tue MaxtMUuTYPE nand 

1h JRELITRNAD EQUATIANS wtTa GECHECT TQ Tet Tv¥Re weasS, TYPE by? 
tus, , TYPE G06 
lAjetae, VVPEnOn9 
DO 1 le1.¥v TYPE N41y 
18s] TYPE 4a 

TV EMBEMOCOVINGT A IRL TODEC Cfo lL Lee 1 Beg, ILPotSPoAVe THe IPIOUe dodo tl oLTVPENG12 

tse TR el)? TVG 1423 
RET YVAN TYhEna14 
ENDC Lode 2s) 1VPE yan5 
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Briss pay 
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7 ee eee el Oe 


Steeit,l 


vive 


« 
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Tr, myth. 


Pe ee a a 
“NT MESV SIP EOD, C1 2%—) BDAC AH, 12Md WEG TYP 
SyeSHEhd WAS Or, 


pe TyPe 
“verb IGE p hal DeVoe UVES, He PS eCOVINCRT YOR 


reettey oot t¥be 


(rn ee a oe oe ee Or Ve? Wl i hk ee 2 er Vewr CORN SET YPE 
MM ML UN ae Oh Hy TAT ADD EL We DPR ee ot Ce ROM, MAT bh T VME 


“MM, 


We Tl pL Thy bebe a Merny ty OB, 


nfot. dpe Cr xs ay, 
ST. US a Vem wt ee 
Ce ae a ho 

~ 1p fob Twbb Ane 

La elhe: 

[Th etme 

mtar fe of Te: 12 Se 
fa gt wye fetta, 
otme deg theyre 

me Tey 

PRG ote «@ wD 


ye, TO ,vibwe,l 


TvOE 


Vet As thr esd TYPE 


AP Ab des 


BPRS OTs 
Aa ACh &,? 


syuersny ¢ 


“Ae tear ade 


Sleek HPT ia em TY Ti or tb teh MATT VRE 


WSBT TO Tae mL MNEs 16 Tek TS Ube ST YP 
aa fees Tyre 
Tybee 
TyPe 


Boom eB btwvl Late te dee CLM Leave OMe 
re oe ee ee ee 
TYPE 
Tver 


Tybe- 
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7ase 
9417 


4:8 
a:9 
a): 
423 
472 
ars 
4/4 
433 


bey ara) 


FuNctlon NEAVETT tue lSelAn i Be TO? TVPE rasa 
WEMerSIO® ved, eo 7CAD, PERE» PERSE7) AVES Id sPIVES Se 7 COVINGSI PE. O3t 
BT bt ERED METESYC 7D GEV, EME oD INC4D OE 120, MD O01 2009 7Od ACT VME, OSE 
DAA AVNET PAY READING DD eSPCOT MOAR ESe I. DOMATEY 062? Tv@Grat, 
COMMON ¥oe7 PEM PEAS, AVEC AVEC: VIND AL BMA, Je TOR eG shots Vo QECORNs SDIYPEr age 
LOM RL IR MTR oD Ryo LM KL RN ENR IAT OU IV, LPENS LU KOAMe MAT SET YPE 635 


APH THOT Ne LT TER, PRAR. METH, Time TYPE @4¢ 
EOUTVALE* Ce ORSAM, vd, Car VATS YD TYPE 443 
THLS ROUTINE COMATES Tre DEVTVATIVe OF Tee [US MOMONT WIth TYPE 445 
RESO-CT TN THE AND CAWAWE TER, TYPE NOS! 
16 Ct A=99909,2 TYPE 10464 
DEAY B mt Tole te te lSote? : TVG 408 
GN TT! 4 TYPE 1) 64€ 
Uc tHe yds, t,6 TVPt 447 
DENVEPERS( | PPERMEL IL, Too tao lS, th) TYPE N44 
IF Ctr t898, oe TYPE 409 
CEAVEDEVeEMELT Lact to TP? TYPE 1450 
6 TO 4 TYPE 1453 
DERVeoPENSE TPDESIGCIIO|N .fao tA ISo IP) TVPE 5452 
LECTPOTSIA, 708 TYPE HESS 
DERVSDERVORIGEL Tuo tao lA. 121 Pd TYPE 1456 
TF 24018d95,0,4% TYPENGSS 
DENyen, Sefeay TVPE 456 
RETURN TYPE 457 
END TvPR nese 
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we Seren 
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1: TEMPOTEMP. CLOVE Pos ede oR dedaVe fal AIEAY (SOL MDI RDERVE 154 oR MoNot Dd 


1? 
a8 


14 BO %) 696, 95015) 5 MER 
9% TEMPSTEMPeANE(T  MdRUUSe(s, 


10 TEMPS TEMPCONE Ce 1 BONER! MPe ht. = ONE Ce LP POCONE CL MIBAVE VoL OKIE 
ONE Ce MdPOAVE Tot, md? 


4? 
qa 


SUQROUTINE NEWTON 


TYPE 1459 


DIMENSION WOS,8r 79,705), RE RE7) PERSEDV oAVE5Ss 7) oCOVES, 507), COVINGSTYPE 9460 
1S eT, AL PwAt DP METERMETIG C7, IRCA IMC OD OVC L2H, 28) 0A 126,126) .00 TYPE Obs 
226), HV27 267, BEPQANC 12D SACSD NOAME De iH rd eMAPE grist, DP 

COMMON X57, PPR, PERS LAV eCOVECOVIN, Al PUA, DE TRAMs GeV shot, Vo BECORN. SDTYPE 1063 
LOMM MKD TA MTR LQs of UNM KL MONO UNK» JATMK CONV, LTEN@ ol Mo MCAMs mal cE TYPE 1494 
APU, IDs STEM, HROWMETH, IMME 

EQUIVALENCE €x©QM, 4), 0G,4ATOY)D 

TIS ROUTINE SETS JP Tre NEWTORCRAPMSON MATRIX EQUATION FAR Tee CHTYPE 1467 
TANGES IN TwE PARAMETERS FOR THE NEXT TTERATIAN. 


MEK eX ARSE (mE TH) 


METH & of AR 9+ [MPLIES SIMPL IE LED Max] MUMeLIMEL THO. EQUATIONS 
METH © ©2 AM oD IMPLIES POMPL ETE MAaNlMUMeLIKELIMN0) FQUATIONS 
METHWO o3 NM oS [MELIES CAMPLETE EQUATIONS “IVIDED BY 9S 


16 (mr Tw) 49,467,8 
Tave. 

DO 27 82,11 
UNSsUC Lede do tee id 
DO 97 yel,ux? 
USSOulLoletods tel? 
00 77 IeJ.ux} 
UTSsulLore tod o€oid 
JaTsi 

TaTalatey 

DO 27 Le2,1Rt 

DO 27 Neiouxi 

DO 27 menN,uni 
JOTOUATo§ 


TEMPODERVI Tyree MeNoL)d 


WCle1) 27,237.90 


GO 7! C1209. 01 omen 


TEMPOTEMPoaViteos MIPNERV Eo 1 om MeMoh D 


16 ¢€ 90194201901" 


VOAV (eg MDAMERV ET, Soe Mo Mot D 


1G CmOLI27 Ld Sou? 
16 (woid16,9601? 


CONE Curt PDO CUNG CUO MdIETSeU0Ss 


LEONE CT MPmAvd ot RPONNEL HePPEAUE Tor em)? 


GO TU 97 


GO Tu 27 
UF 6901927,97018 


TEMaCOVE at Mor eLOECAVE vel Ned LIOHOUEC TO TereLel DaCNu( Jyorzemegoey) 
GN TO 62 6 1O01Od MER 


TEmelEmsuneg 

16 (MOND 22, 91022 
Téme: .SelGu 
TEmpaTENPoTEM 
ACLAT VAT IATEMP 
Le 

00 38 Kes, 18) 
Letei 


TYPE 1462 


TYPE 405 
TYPE 9466 


TYPE y468 
TYPE 0469 
TYPE ne79 
TYPE O74 
TYPE N472 
TYPE 94673 
TYPE 9474 
TYPE 5475 
TYPE 0476 
TYPE NG77 
Y¥0G 9478 
TYPE NG79 
TYPE 480 
TYPE 5604 
TYPE 402 
TYPE 9483 
TYPE (604 
TYPE 9405 
TvPE 406 
TYPE N48? 
TVOE n4e8 
ALT OEY } | 
TYPE 9490 
TYPE NGOs 
TvVPE 402 
TYPE 5403 
TVPG 404 
TYPE 0405S 
TYPE 0496 
TYPE 3407 
Ade DOT | 
TYPE 9609 
TYPROSCO 
TVPG 081 
TYPE NS08 
TYPE 0903 
TYPE S04 
TYPE 0505 
TYPE 0906 
TYPE 0807 
TVPE 0908 
TYPE 9909 
TYPE 9S10 
Ald Tides 
TYPE 5912 
TYPE S15 


Bf BIL ISUNSE COV ToL UeL eK eAVE los oMIAAVE SOL oR DOAVES o3.Kd8 
WE pologqodomnLIOAVe JoqeRdmiil sng, 10) ohogd 
33 


31 
38 


34 


UOSsUCL.101510€01) 
BtLrver1. rs 

DO 29 182,mx4 
Letel 
BcLyeavile,m) 


“GO Th €29,98229),.MEK 


OcLreBeL res 

BEL v@Bch denote Le toKed? 
Do 33 ye2.ux1 

DO 33 leJ.wyi 

LeLet 

BO TU 632.3: 63: >. mMEK 


60 T0 

BELIOCOVE Pate deLeMIOAVElLa~1 KI MAVE Jol eK) 

StL peWch dere tele doKod? 

GO VO 637437034). MEK 

faven 

DO 36 Ke2,1R1 

JATou 

DO $5 Le2,1A1 

DO 35 Neot,wx1 

DO 35 MeN, mx 

JATaJATo2L 

OVC VATIS DOAVELe Le MoMaNoLIAUCLo Lele loModd 

OO 3@ Jei,uni 

OO $@ fed, uxt 

SaTalates 

JATe0 

DO 34 Led, thi 

DO 36 Nei,mys 

00 36 MeN, mNd 

JATaJALoL 

ACLAT SATIS ALCTAT JATIOBECLATIOBVCUATS 

RETURN 

00 44 384,98 

DO 43 vei.t 

ACL VIOUCLedodo le toinvet? 

ACJ deatt we 

Ole UlLoeatedotodtelosres.o 

CALL MATINvEA,1850.2.08) 

LeieR 

OG 92 mod, 1R1 

PERTS UlLegedstedetd 

BO $2 lei,uy 

LeLed 

SCL IOUC Lot ete lode re RISPERT oAve TK) 

Ovitde AVET RdeReL) 

BO $2 Jei,t 

LOLo4 

Steen tedene tyke ne Loh {Pent OOVELDeRY Cd COVE Te deK) 
3 


@0 fn 
wD 


TYPE 5534 
TYPE 4595 
TYPE S36 
TYPE 0837 
TYPE 596 
TYPE 1545 
TYPE DS2r 
TVPE 9821 
TYPE 9$2% 
TYPE QS22 
TYPE 9924 
TYPE $28 
TYPE nS2e 
TY°6 9927 
TYPE nS 20 
TYPE pS29 
TYPE 033¢ 
TYPE 0933 
TYPE OSSE 


ie 


-~@a~w 


ry 
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9” 
191 
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ee 
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Sus? UTINE RAPwSON TYPE 9968 
DIMENSION 065,30::790266), PERC) PERSE 7) DAVES 57) oCOVIS, 507) COVINISTYPE 1 9G9 
LoDo PP. ALPwAl 7)  DETERME7) GOT, TKE4Do TNE420V0 122 289006 126,190) Be TYPE OS 74 
226), 840124) RECNANI12)- SACS) MEAP (Sogn oMalhinneti ? T¥PEDS72 
COMMON KL7, PER, PERS, AVO COV. COVIN, ALPHA, DETERM,G.VoAod, RVs RECORD s SDT VPE (S72 
10ME RG. 18, Ia. Tq. LUMNE MLR OND RNKs JATMN CONV, TTERM LK Keam. eat oETYPE (S79 


2Pm,TNe ID 1 TER, PROD, METH, 1 DLIMP TYPE n974 
COULVALENCE (xAmM, vd, Cae maToN} TVPE 9375 
Veale ROUTINE ADDS THE 2PPROPAL ATE CmangES TO THE PakamETERS To OCOTTVPE D976 
JAIN PHBIR VaLURS FOR Tue NEXT 1TERATION. TVPG NS?) 
anare?,: TvPE 970 
UF cme MIO, 30.23 TVPE 9979 
WEWYONCRAPUSON IT TERAATION TVPG 090C 
saves TYPE 9Ay 
OnTotmugecunyesd)72 TVPG 99R2 
SwoRVEn INCREMENT VECTA® UNTIL ALL PENCRNTAGES ARE WITHIN BOUNDSTYPE (S03 
00 3c we2.1R: TvOE 0904 
AMAPOSPERS (MISC 2. OC TAT» TVPG 9909 
16 camayy7. 7," : TvPEn9ee 
QAP CAKb7 © 8 SsActat? TVPG 0907 
UF Cam apodwaldO.2 3.49 TYPE 0900 
Amaveawaz TYPE 9909 
lavelavomyy TYPE 990u 
19 (am hog .nd21.13.13 Abe Ti Le | 
BO 22 ned, Aran TySG09C8 
Oemreamavaacn) TYPE 09O3 
tat ef TYPE 0904 
80 £0 me2,182 TvOG 0908 
tatelates TVR 996 
PERRICK OPENSIM OAC IAT) Tv0G 4907 
BO a4lort my Ty¥PE SCO 
latealares TVPR 0900 
AVEToMP@AY eT oMVOMETAT) = TYPE 0000 
00 99 Jes,ux 1706 0603 
DO aSledsMy TvPb0002 
tatelater TvPR 0003 
COVINGL ce MPOCAVING os QoBtomtiard TVPENO04 
CAVING Se DoW POCAVING Ts Uek) Ty¥®G 009 
BCE SPOCUVING Te JoRD tyPBne0e 
Otsoldeatt,v) Tvb1607 
CALs MOVinvc dwn Ay.* as TVPG 0008 
AeahsledaysDa TVPR G09 
DEH MCR Pe CRONTE Calenarread TvPE O10 
10 ¢01997027.28 TvPErea1 
Inve rene, TVPG 0612 
URTTE NUTEINT TAPE 25171, TNLEY, NA - TvOB 619 
CHMRAT (Ade HETEMMINANT CE COV ADTANCE INVERSE FG Trem Tes mek ya, Ort vVPE res ¢@ 
DETEMMITDe o:. TvPG O13 
DO 2@ yei.ur TvPb O16 
DO 96 lel.un TvPbne17 
COVE lsd dealt ou? TVPG 618 
RETION TYPE 619 


toy 


$3 


SUCCESSIVE SUBSTIT TIONS, 


- LaeyR 


Mx a) 


SMORTEN INCREMENT VECTOR UP TeL ALL PERCeTAGES ARE wWTTMrh BOUNDS 


DO 34 «87,181 
AKA7E<PERSCKISC2,¢ SRK TDD 
TF CAMAPIST, 33, 36 

AKA? & AKAZ ©: .S5S/A( KO) 

IF Cam bPodeaty 3S, 34,34 
AMATOAKAZ 


hb CONTINUE 


LF CAKAT C3 ond87,39, 39 

Bo 38 wey.1A 

UCKyeARAT UAH) 

DO 58 We 4, 1h) 

PERS CK ISPERSIKRDS ALKod) 
DO S2 3ei,mx 

LeLe.t 

AVOT SRD SAVE] oRVCH(L) 

DO 82 Jei,! 

Lele 

COVED DoKMCOVEL se Ved OMILD 
COVE sol eM dRCOVEL 0K) 

ACE JieCOVEl» U.K? 

ACJ, TdeACT, 0? 

Con? thus 

CALL MATINved. Mx, Ov, 08) 
ApsadSF (Dar sna 
DETERMCKDSCRARTECANDADDeAD 
DO 84 tes,ux 

BDO 84 Jyel.uH 
COVING Dede WIOATT ou? 


a CONTINUE 


60 TO 29 
END 


a7 


TYPE 620 
TYPRnAZL 
TYPE 3622 
TYPE 1623 
TYPE 1624 
TYPE nO25 
TYPE N626 
TYPE 5627 
TvP6 4698 
TYPE 629 
TYPE 0636 
TYPENO633 
TYPE nO32 
TYPENOSS 
TYPE 5034 
TYPE 695 
TYPE NO 36 
TYPE 0637 
TYPE NOSO 
TYPE 639 
VVPEnO40 
TVPG 0643 
TVPE 6468 
TYPE 043 
TYPR read 
TYPE 0045 
TVPENO46 
TVPE 0647 
TYPE nO48 
TV¥PE0O49 
TVPE06SO 
T¥PE0OS1 
TYPE0SS2 
TYPE DOSS 
TYPE 9694 


ba 


wT) 


SUBRUUTINE RESULT TYPE n655 
DIMENSION ¥05, 30 G702716), DERE 7) PERSE 7) AVS. 7) sCOVI9, 5.7) -COVINCSTYPE 1656 
12 Se Fo ALPMALT)  NETEMMETD,REPV IRC OdO INC 40 V6 12% 7890861 26,1260) Bt TYPE 9657 
P2EdHVEZ IAD RELNAN( 32) 2 SNCS) KCAMESoi gn oMaT dorset, ? TYPE 1658 
COMMUN 27, PER. PERS, AV OCOV,COVIN, ALPHA, DETERM, Ge Voie 4onV oRECORNs SDT YPE 1659 
LoMM MN TRIMI Re Lae | UMNK LAL ONY NK o JATMM CONV, J TERMoL YY -MOAMe maT ETYPE n66n 


2PML EXO IN TTER,PROR, METH, I NUME TYPE n663 
EOUTVALERCE EXSAM,V),¢A, MATA) TYPE (062 
TWES ROLTINE PRINTS QUT Twe SCAQAMETE®S VESCRIBING EAC TYDE. TYPE 1063 

1 FORMAT C 44, 48m, SAMMAV IML MOL IME! TmOOD ANALYSIG NF TYPESs/Gv,9346) [YPE 664 
2 FORMAT C11, 1 5R > SOMCMARACTERISTICS NF THE weOLE SAMPLES //3u) TYPE 5665 
BS FORMATCSNSS/A// 40> 2 SHCmMARACTERISTICS AF TvPE 14as73x) TYPE noe*% 
@ FORMATOLIN, 49KTHE PROPORTION AF TwE POPULATION FROM THIS TYPE of. 1965667 
1373) , TYPE cee 

5 FORMATOCS. xX, SMMEANGS/5&6 12.9) TYPE coe9 
& FORMATERESZONN, 1OMETANDARN NE TATIONS s5F 12.2) TYPE 0670 
7 FORMATO URJ9AN, L2MCORRELATIONS) TYPE 0071 
a Format (5649.4) TYPE O72 
@ FARMATOSNASAAA 120, LSHSAMB LD SIFE SIL IS LUMO PIMNUMHER OF VARTAOLES ef YPE 9073 
D197 ML LP MBER AF TYPES s14/ 725K, 1AM TTERATION NUMMER [3//39K, TYPE (074 
PISMLTRELTMAND IF 12,230 TYPES FN THIS SAMPLE 9648.8) TYPE 0678 
0O 3" veleux TYPE 0676 
Jasyel TVPE O77 

Vt AVG Se LIMUC Le Lele Jdore*? TYPE ne78 
DO 2) Jaloux TYPE 0679 
Jase} TVP6 0600 

DO 2). Lade¥x TYPE 602 
laaje. TYPEn6e2 
COVED eo Teld aire pe tbo tho let IobVedeod SAVES o1? TYP&0OGS 

9 COVE DeLee TAVET ede? TYPE 00A4 
WALT NUTHUT TAPE 2,7, ¢RECORN ETD, 1830949 TvPE 085 
WRITE NUTHHT TAPE 2.ONK, MHD IR, LITER IR, PmUM TVPE 1686 
QO +h wet,1a) TYPE DOA? 

Ki sKe} TvPé 088 

1F cad) 49,912012 5 TYPE ..609 

41 WRITe OUTPnt TAPE 2,2 TYPE .,090 
Go Tr 43 TVPENO02 

12 WRITE QUTPHT TAPE 2.39.K1 TYPE 602 
WRITE AUTHIT TAPE 2,4, PERSERK) TYPE O93 

qt WRITE HUTeNy TAPR 2.5,CAVCT Od Tere me) TYPE 604 
00 14 Jeieux TvPE NGOS 
SDveCOvel.t,") TvPE 0006 
AfeanMFeanyr seny TYPE 0607 
SAve(SAATF aNeeovedean ; TvPE Eee 
o@ SNCTPsSpV TYPE 0609 
Of 19 Teieun TYPB 0700 

GN 1SJeloMy TYPE 703 
ACE DPOCUVE Te reRDFESNEL DASE C DD TYPE 0702 

a8 ACS Tieatys TYPE, 703 
BRITE GUIPHT TAPE 220,05 CT De Teiemnd TVPE5 704 
URIT NUTOnT Tape 2,7 TYPEN7NS 

UN 1A Jateuy TYPEN7nO 

qe WRITE OUTP TE TAPE SoS CAEL esd, 1Sqomud TYPE) 7097 
RAETimN PvP 0708 
EN ols ore.) TVPR (708 


c 


SUARUUTIFE PILACe TYPe0710 
DIME .STO® vet5,%- de 7LOD, PERE IV, PERSE) o AVES. 7). COVIS, Ne 7d oCOVINESTVPES 711 
tS Po AL PWAC TY, METERME ID GEIY, EMCO DIMES) 102 2h5 2800 M11 20,120) OCT TYPEnI 12 
DOAD, AVE Phd, MEPNAMC DDOSMESD UCAMCS rte DoMETEY 61. 9 T¥PEN713 
COMMUN O27, PER PERS AVeCNV eC OV ING ALPHA Ohl sGov ode t eb oe RECORNs SDT VPEL 714 
TMM, AMT LR MT Re TR Te LIMNE KE ORO NM ATER s CONV G TTERY of Me ROAM eo MATS ETYPR NZIS 


PPuwet eo INe TTER, CROR MET HAITI ME TVPE 4746 
FQUIVALESCE (XSAN, vd, Cae MaTOYD TVPG N77 
THIS ROUTICE PRINTS AUT THE CRABOBILITIES OF TYPE MEMAERSUIP FAR ETVPE 1778 

TACH OBSERVATION TvPE0719 
FORMATCIN) 9%, TOHPROQANTE TITIES OF TYPE MEMNENSHIDS//4, 7184 TYPE N720 
FORMATO} ¥, 764, %) TvPEH723 
WRITE OUTPVT TaPe coi, Cloteredwd TYPEN722 
KOs" : ¥VPE 723 
SRENSD TYPE Q724 
REwtD 4 TyvPE (725 
CALL BLOCK EAH, 81 Hh KRMNWSTSRERD SL ONG) TYPE 0726 
REAM TAPE a, CONCTe do Lar eMED, 18101 ONG) TVPEQ727 
DO A Kayel AVG TYPE N7 26 
DN 4 leademy TVPE:: 729 
Zelotdewet,«? T¥PE Y73H 
Catt DENSITY TYPE O733 
KARanhkes vPEN732 
DO ® 1ler,tes TvPEN733 
GelryeGcldeeERsel) TYPE N734 
WRITE HUTT Tame 2.9-KO, CGC 1e2.1A9) TVPEN7SS 
16 ¢609709,3 T¥PE 1736 
RETUAN T¥PG 737 
ENDE-oLetionor? TYPEN736 


AAA 


SUBROUTINE CCK IND LALK, KRMN, AST OKEND LG WD 

TmIS ROUTINE COMPO TES ViwgeRS iSeFUL IN Tre CONTROL OF Tee TNBUT 
AND OUPTRUT OF LISTS MF LENGT® WX IN ALACKS OF LEAGT™ LULK, 
KRANONUMHER OF [TEYS VErPAINING . Tee Last 

LONGOLENGTw OF CURSERT BLOCK 

KST AND KEND ARE TRE STAQTING AND FNOING IADEMES FOM Tue FTEMS IN 
TwE CURRENT HWLACK, 

1F (RRENI2Z, 102 

KRMNONX 

MENDED 

LONGSLALK 

KSTeKEN/De1 

TE CRMMNOLAL OPT, 406 

LONGSKRMN 

RENNOKENL oO: ONG 

HAMNOKK MN], ONG 

RETURN 

tnd 


TYPE: 739 
TYPG i74au 
TvP@: 741 
TVPG 742 
TYPE D743 
TVPG17446 
TYPE 3745 
TPR 1746 
Ty@e3747 
TYPE! 748 
TyPe749 
T¥eeu7Sh 
TYPEu7S1 
TyPk 752 
TYPE 753 
TYPE 754 
TvPe..755 
TYPE 796 


AAAaEAAANUAANAAAN 


SUAROUTINGE MATINVEALN, Bow. LETERM) 


DIMENSION TOTVATC1 26), R01 2001200 BEL 2Oet Ie INNEN(1 2602) PIVOTC126. 
EQUIVALENCE ¢(180W, SRW). CICOL IM, VCO UMD. 


camax, T, 


SwaP) 


PROGRAMMED HY BUNTON S. GaRBOW, ARGOA"E NATIONAL LABORATORY, 
AND PEPOSTED Iv [ay 7 627:9 SHARE LIWRATY aS AN F4 ¢, 

TwiS SUARNUTINE CIMELTES Tee [NVERSE 4NN DETERMINANT OF 
MATRIX As AF CADER A, HY TRE GaUSS-UARDAN “EPHOD. AcINVERSE 
REPL ACES A. AAD THE NETERMINANT CF 4 TS PLACED IN DETERM. IF 
Ms 1 THE VECTOR @ CONTAINS tre CONSTANT veCTOR wHev MaTinv IS 
CALLED, ath TRIS TS REPLACED witm Tee SOLUTION VECTOR, IF me 4, 
NO SIMULTANEGOLS EQUATION SOLUTIONS AME CALLED FORs oND BW PS NOT 


PERTINENT, N IS NOT TO EvCEEn J. 


As Ne Wo M. AND VETEQM [uy Toe gRGUMENT LIST ae PUM™Y VARTARLES. 


INITIAL IZA ON 

DETEAMe1.0 

DO 2uv Jsi.® 

iPLvOTCude4 

DO $50 te1,N 

SEARCH FCR PIVAT ELEMENT 
Amayet,: 

00 105 yei.N 

16 CIPTVOTE 991) Bue 108, 69 

DO 100 Ke1,N 

1€ CLPPvVOTemde1d Bue Tei, 74" 
16 CABSE CAMA ADSE CAC Te WddD Bus 
TROueu 

(COLUMeK 

AMPAWOAC JON) 

Cont! Nué 

CONTINUE 

UPTVOTEICOL UMPETPTVOTC ICAL UMS os 


love 


210 


NTEACHANGe ROWS TO PUT OFVOT ELEMENT OW DIAGONAL 


1F CEROweT ALUM) 1465 26 » 149 
DETEANeoDETERH 

00 200 Le2,N 

SuaPoAcTPOuet? 
AC1ROwWee dear ICOLumet) 
ACICOLUM.t pe8uwaP 


TYPES 757 
TYPE" 758 
TYPE, 759 
TYPE: 7a0 
T¢Pe 702 
TYPE 742 
TYP 17635 
TVPG 704 
TYPen 765 
TYPeL 766 
TYPE0767 
TYPE 708 
TYPE1709 
TYPEr 7 IC 
T¥OG 771 
TVPE 772 
TYPE 773 
TYPE. 774 
TVPE5 775 
TVPE 770 
TY@G)777 
TYPE 778 
TYPG0779 
TYPE 70U 
TYPGu 762 
TYPO. 702 
TVPEU 703 
. OEr 764 
TYPE" 765 
*7PG0 786 
Tv 53767 
TYG ° 788 
Tree. 769 
TVPE”, 90 
TYPENHIV3 
TYPED 72 
TYPEG 703 
TYPE. 704 
Tv0G5 708 
TVPEC 796 
YvPG0707 


c 


TF Cmd 26r, 261, 22. 

OM 230 L981, ™ 
SviaPoOc1@Ou..? 

Oc ROW.) PemcICoLum.L?> 

Oc lCOLumst »eSusP 
INDER(C I. 2201 R80w 

twDEKil. 29<!COLuM 
PIVOTCEDeAeTCOLUM. 1COL UM) 
DETENMeNE TERM aPIYATIT) 


DIVIDE PIVAT ROW BY E'VCY ELEMENT 


ACICOLUM, FeAl; 04.0 


1 DO 390 Ler 


ACICOLUMeL POALCICOL UML IsO1.OTCT) 
(Kem) SBr. BA:, JA. 


» DO 37m Leor,m 


Bc lEoLumst veBECOLUs..dsPTVOTCEL. 
REDUCE NONePIVAT BONS 

DO 5$c Lie1.N 

Fer respCMeemds 42), $8, @ . 
Teathi,1ca um) 

ACLaclCOLiwsre 4 

DO 496 Let.n 


ib ACL oL DOA *eLPCACICALUM, LIB? 


16cm) $a:, as e 4A. 

DO 5:5 Let,m 

BCL oL deb seLVoACICOLUM,L Ie? 
CONTINUE 

ImTENCHANGE COLUYNS 

OO 720 Tet 

LeNneie} 

16 CUNNGr es LOC TNMEMEL 29) Of, 
JROWSINDE NCL 21) 


* JCOr Mat nnaxtt. 2? 


00 7.5 Ket,N 


» SUAPSAIR, la0e) 
ACK, IMOwI SOE. UCM Um) 
* ACH, JOM Ure Saar 


CONT [ee 

CONT [nue 

WG Tisen 

GND 

ENNG adele te D 


7L4: 0 O3y 


TYPE 0708 
TYPE R709 
TYPE DRC 
TVPG 0OH1 
TYPE QOn2 
TYPE 9603 
TYPE 0004 
TYPE 9809 
TYPE 4896 
TVPE 0097 
TYPE 9808 
TYPE 9009 
TYPE 010 
Mla rite! 
TVPG 0012 
TVPG 0813 
TYPE 0014 
TvPG 0019 
TYPE 0036 
Tyv?G 0037 
TVPG 0018 
TYPE 0019 
Tv9G 0020 
TvP6 0021 
TvPE O22 
TYPE OO23 
TYPE nOa4 
Tv0G 9098 
Ala ti ta] 
TvPG0G2? 
tv06 0020 
tveG eed 
TvPGHOS0 
TVPG O31 
TVPG 0632 
TVPG 039 
TVvPG 0034 
tveg ness 
Tv0S 1036 
¥v000037 
TYPE yO38 
TYPE NG SO 


